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A New Method for Recognizing Quadric Surfaces from Range Data and 
Its Application to Telerobotics and Automation 

(phase II) 

by 

Nicolas Alvertos* and Ivan D’Cunha** 

Abstract 

Pose and orientation of an object is one of the central issues in 3-D recognition 
problems. Most of today’s available techniques require considerable pre-processing, 
such as detecting edges or joints, fitting curves or surfaces to segment images and try- 
ing to extract higher order features from the input images. In this researchNve present 
a method based on analytical geometry, whereby all the rotation parameters of any 
quadric surface are determined and subsequently eliminated. This procedure is itera- 
tive in nature and has been found to converge to the desired results in as few as three 
iterations. The approach enables us to position the quadric surface in a desired coordi- 
nate system, then^ utilize the presented shape information to explicitly represent and 
recognize the 3-D surface. Experiments were conducted with simulated data for objects 
such as hyperboloid of one and two sheets, elliptic and hyperbolic paraboloid, elliptic 
and hyperbolic cylinders, ellipsoids, and quadric cones. Real data of quadric cones 
and cylinders were also utilized. Both of theses sets yielded excellent results. 


* Assistant Professor, Department of Electrical and Computer Engineering, 

Old Dominion University, Norfolk, Virginia 23529-0246. 

** Graduate Research Assistant, Department of Electrical and Computer Engineering, 
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1. INTRODUCTION 


One of the most important tasks in computer vision is that of 3-D object recogni- 
tion. Success has been limited to the recognition of symmetric objects and now 
researchers are investigating recognition of several asymmetric objects as well objects 
placed in complex scenes. Unlike the recognition procedure developed for intensity 
based images, the recent upsurge of several active and passive sensors extracting qual- 
ity range information, has lead to the involvement of explicit geometric shapes of the 
objects for the recognition schemes [1]. Location and description of 3-D objects from 
natural light images are often difficult to obtain. Range images on the other hand give 
a more detailed and direct geometric description of the shape of the 3-D object. 

The determination of the location and orientation of a 3-D object is one of the 
central problems in computer vision applications. It is observed that most of the 
methods and techniques which look into this problem require considerable pre- 
processing such as detecting edges or junctions, fitting curves or surfaces to segmented 
images and computing high order features from the input images. Since 3-D object 
recognition depends not only on the shape of the object but also the pose and orienta- 
tion of the object as well, any definite information about the object’s orientation will 
aid in selecting the right features for the recognition process. 

In this research we put forward a method based on analytical geometry whereby 
all the rotation parameters of any object placed in any orientation in space are deter- 
mined and eliminated systematically. With this approach we are in a position to place 
the 3-D object in a desired stable position thereby eliminating the orientation problem 
and subsequently utilize the shape information to explicitly represent the 3-D surface. 

In the initial part of the research we discuss the rotation transformations and in 
the later part of the research we propose our scheme to eliminate the rotation 


parameters. 


2, BACKGROUND 

Most of the available techniques for describing and recognizing 3-D objects are 
based on the principle of segmentation. Segmentation is the process in which range 
data [1] is divided into smaller regions (mostly squares). These small regions are then 
approximated by planar surfaces or curved surfaces based upon the surface mean and 
gaussian curvatures. Regions sharing similar curvatures are subsequently merged, the 
process known as region growing. There are several other approaches [2,3,4, 5, 6] 
wherein the 3-D recognition problem has been dealt with. Levine et al. [7] briefly 
review the various works in the field of segmentation, whereby segmentation has been 
classified as region-based and edge-based approaches. Again surface curvatures play 
an important role while characterizing each of these approaches. 

We have proposed an approach [8] based on two-dimensional analytic geometry 
to recognize a series of three-dimensional objects. An effective technique to determine 
the 3-D object location and orientation will aid us in extending the proposed scheme to 
various 3-D objects such as the hyperboloids of one and two sheets, the circular and 
elliptical quadric cones, the circular and elliptical cylinder, the parabolic and hyper- 
bolic cylinders, the elliptic and hyperbolic paraboloids and so on. Figure (1) illustrates 
the various 3-D surfaces we propose to use for the recognition scheme. In our pro- 
posed method [8] a feature vector consisting of various 2-D curves obtained after inter- 
section of objects with planes in various orientations serve as the medium of distin- 
guishing objects from one another. Sets of range data of objects are attempted to be 
identified as quadric surfaces based upon their representation by a second degree poly- 
nomial. 




Figure 1. Quadric surfa 
cone. Hyperboloid of or 
bolic cylinder, Hyperbol 






Recently [9] recognition of 3-D objects based upon their representation as linear 
combination of 2-D images has been looked into. Transformations such as rotation 
and translation has been considered for the 3-D objects in terms of the linear combina- 
tion of a series of 2-D views of the object. Intersections of various 3-D objects with 
each other, the edges obtained thereby, has been [10] looked into as means of obtain- 
ing a surface description vector (SDV) graph for representing natural quadric objects. 
This research proposes an optimal sensing strategy of a sensor system engaged in the 
recognition and localization of 3-D quadric objects. 

3. THEORY 

Any quadric surface can be represented in terms of a second degree polynomial 
of variables x, y, and z, such that 

F(x,y,z) = ax 2 + by 2 + cz 2 + 2fyz + 2gzx + 2hxy + 2px + 2qy + 2rz + d = 0 (1) 

Let (x,y,z) describe the coordinates of any point in our coordinate system. Con- 
sider a rotation of a about z axis, then the old coordinates in terms of the new one 
are represented as 


x = x'cosa + y'sina 

(2) 

y = -x'sina + y'cosa 

(3) 


i.e., the rotation matrix is 


cosa sina 0 


Ra = 


-sina cosa 0 
0 0 1 . 


(4) 


Next consider a rotation about the x' axis by an angle p of the same point. The 
resultant coordinates and the old coordinates are now related by the following equa- 


tions 


z = z'cos(i - y"sinP 

(5) 

y' = z'sinp + y"cosp 

(6) 


whereby the rotation matrix is 

'l 0 o' 

Rp = 0 cosP sinP (7) 

0 -sinP cosP_ 

And finally consider a rotation about the y' axis by an angle y, then 


= z"cosy + x"siny 

(8) 

= z"siny - x"cosy 

(9) 


The rotation matrix for the above process was 

cosy 0 siiry 
R Y = 0 10 

-siny 0 cosy 

Figure (2) illustrates the behavior of a particular point (x,y,z) with respect the various 
rotations described above. 

Using the above set of equations and solving for x, y, and z in terms of x", y", 
and z" yields 

x = — x"(cosacosy + sinasinpsiny) + y"sinacosP + z"(siny:osa + cosysinasinP) 
y = x"(cosysina - sinysinpcosa) + y"cosPcosa + z"(cosysinPcosa - sinysina) 
z = -x"siny:osP - y"sinp + z"cosycosP 


X 



Figures 2(a), (b), (c), and (d) refer to the coordinate system initially, 
after the first rotation, the second rotation, and subsequently the third 
rotation, respectively. 



After substituting these new x, y, and z coordinates in equation (1), we get an entire 
set of new coefficients for x 2 , y 2 , z 2 , yz, xz, xy, x, y, and z. These new coefficients 
are as listed below : 


A' = cos 2 / [a-cos 2 a + b-sin 2 a] + sin^sin 2 / [a-sin 2 a + b-cos 2 a]+ 2sinasinPsinycosacosy (a - b) (10) 
+ c sin 2 )cos 2 P - sin2a [h sin^sin 2 /] - sin2y [fsinacosp + g-cosacosp + h-cos 2 asinp - h-sinPsin 2 a ] 

+ sin2Psin 2 y (f-cosa - g-sina) + h-sin2acos 2 / 

B' = (a-sin 2 a + b-cos 2 a)cos 2 p + c-sin 2 p + sin2P [-f-cosa - g-sina ] + h-sin2acos 2 p (11) 

C' = sin 2 / (a-cos 2 a + b-sin 2 a) + (a-sin 2 a + b-cos 2 a)cos 2 ysin 2 p + 2sinasinPsinycosacosy (a - b) (12) 
+ c-cos^os^ + sin2a [h-cos 2 ysin 2 p - h-sin 2 /] + cos 2 ysin2p [f-cosa + g-sina] 

+ sin2y [-f-sinacosP + g-cosacosP+ h-cos2asinP] 

2F' = j^(b-cos 2 a + a sin 2 a + h sin2a - c)sin2p + (2g-sina + 2f-cosa)cos2pj cosy 

+ j^((a - b)sin2a + 2h-cos2a)cosP - (2g-cosa - 2f-sina)sinpj siny (13) 

2G' = sin2y[-cos 2 a(a + b-sin 2 P) - sin 2 a(a-sin 2 p + b) - c-cos 2 p (14) 

- sinPcosP(f-cosa + g-sina) + h-sinacosacosP 

+ sin2a(a + b)sinp + 2cosPcos 2 y(f sina - g-cosa) + 2h-sinP(sin 2 asin2y - cos^cos 2 /) 


2H' = sin2a cosacosP(b - a) - h-sinPsinycospj + sin2Psiny(a-sin 2 a - b-cos 2 a + c) 
+ cosysinP(2g-cosa - 2f-sina) + sin 2 Psiny(2g-sina + 2f-cosa) - 2h-cos 2 acosycosP 


2 F = 2cosy [-p cosa + q-sina] - 2sinPsiny [p-sina + q-cosa] - 2r sinycosP 


2Q' = 2cosp [p-sina + q-cosa] - 2r-sinp 


(15) 

(16) 
(17) 


2R' = 2cosysinP [p-sina + q-cosa] + 2siny [p cosa - q-sina] + 2r-cosycosP 


(18) 


D' = D 


As seen from above, except for the constant D, all of the other coefficients are 
affected by the rotation of a, and y. 

3.1. Scheme 


The product terms yz, xz, and xy in f(x,y,z) above, denote the rotation terms 
which are to be eliminated. Elimination of all these rotation terms will place the 3-D 
surface on a coordinate system plane parallel to our coordinate system. 

At first sight, this problem appears quite simple. Eliminating each term, i.e., by 
rotating the surface about the origin in a particular orientation by a suitable angle 0, 
then eliminating the second term and then the next term. But that’s not the case. In 
the presence of a single rotation term i.e., if the equation is in the form 

F(x,y,z) = ax 2 + by 2 + cz 2 + 2fyz + 2px + 2qy + 2rz + d = 0, 

equation of the trace of the above surface in yz plane is obtained by setting x = 0, so 
a rotation about the origin in the yz plane by an angle 0 will eliminate the term yz 
from the above equation. 

However, in the presence of two or more rotation terms, while trying to eliminate 
a second rotation term, it is seen that the rotation term eliminated first, reappears and 
thereby at any given instance there will be at least two rotation terms available. The 
approach we have proposed is an iterative process, whereby at each stage the object is 
rotated in each of the directions x, y, and z sequentially and the process is carried out 
until all the product terms fall by, i.e., the coefficients f, g, and h converge to zero. 
Since our aim is to eliminate the rotation terms xy, yz, and xz, let’s now exclusively 
consider the coefficients of these rotation terms, namely F, G, and H evaluated above. 
In an iterative procedure we will be able to eliminate all of the product terms. For 
e.g., we wish to eliminate the term xy. Then, by a specific rotation of a about the 


z axis, we will be able to accomplish our goal. However, while carrying out this pro- 
cess, the orientation of the object about the two planes yz and xz , i.e., the angles 
the object made with these two planes have been disturbed. Now if we wish to elim- 
inate the yz term, the object has to be rotated about the x axis by an angle p. How- 
ever, in this instance, while carrying out the process, the already missing xy term 
reappears though the magnitude of its present orientation has been reduced. Hence by 
carrying out the above process in an iterative fashion, there comes an instance when 
all the coefficients of the product terms converge to zero. 

Let us consider the equations (13), (14), and (15) respectively. Let us eliminate 
the coefficient h, i.e, the xy term in step 1. This can be accomplished by rotating 
the object about the z axis by an angle a, whereas p = y = 0. Under these cir- 
cumstances the new coefficients look like as shown below: 


2f u = 2g-sinai + 2fcosa 1 
2g n = 2g-cosa 1 - 2f sinai 

3 

2h n = (a - b)sin2a! + 2h-cos2a! = 0 where cot2a! = b - — 

2h 

As seen above the coefficient h has been forced to 0. The most significant bit of the 
subscript refers to the iteration number , whereas the least significant bit of the sub- 
script reflects the number of times the object has been rotated by a specific angle. In 
the above case the LSB of 1 refers to the first instance where the object has been 
rotated by an angle a. The remaining coefficients a, b, c, p, q, and r also reflect 
changes brought about by the above rotation. The new coefficients look as shown 
below: 


a n = a-cos'T*! + b-sii^oq - 211 -sinc^cosai 


b n = b-cos 2 <Xi + asin 2 a, + 2hsina 1 cosa 1 


c„ — c 

2p u = 2p-cosa, - 2q sina] 

2q u = 2psina 1 = 2qcosa, 

2r j j — 2r 

The new quadric equation now has a look as shown below: 

F(x,y,z) = a u x 2 + b n y 2 + Cnz 2 + 2f n yz + 2g u xz + 2p n x + 2q n y + 2r u z + d = 0 

Consider the second step wherein the coefficient corresponding to the yz term is 
forced to zero. In this particular case, the object has to be rotated by an angle p 
about the x axis, whereas a = y = 0. Under these circumstances, the new rotation 
coefficients (signifying the product terms) becomes 

2f 12 = (b [2 - c 12 )sin2(3i + 2f u -cos2p 1 = 0 where cot2Pj = 11 — — — f u 

2g, 2 = 2g u -cosPi 
2hj 2 = -2g n sinPi 

At the same time the other coefficients become 

a 12 : a U 

bi 2 = c u -sin 2 Pi + b n 'COS 2 Pi - 2f n -sinp,cosp 1 


Cj 2 = bifsin^! + c u cos 2 Pi + 2f 11 -sinp 1 cosp 1 


2p ]2 - 2p u 


2q 12 = 2q n cosPi - 2r u sinP, 

2r 12 = 2q n sinP, + 2r n -cosP, 

The new quadric equation as before looks like as shown below: 


F(x,y,z) = a 12 x 2 + b 12 y 2 + c 12 z 2 + 2g 12 xz + 2h 12 xy + 2p )2 x + 2q 12 y + 2r 12 z + d = 0 

In the final step of the initial iteration, the coefficient corresponding to the xz term is 
forced to zero. In this case, the object is to be rotated by an angle y about the y 
axis, whereas a = p = 0. Under these circumstances, the new rotation coefficients 
become 


2f,j = 2h 12 -sinYi = — 2gj [ •sinP 1 siny, 


2gi 3 =(a 13 - c 13 )sin2Y] + (2g n -cosa! - f u sinai)cosP]Cos2Yi = 0 


where 


cot2Yi = 


c 12 ~ a 12 
2 


Sl2 


2h 13 = 2h 12 cosYi = -2g 11 -sinP 1 cosYi 

Let’s now carefully analyze the coefficients of xy, yz, and zx obtained in the final 
step of the first iteration. Consider for instance the coefficient corresponding to the yz 
term. It is observed that while proceeding from one step to the other, the new 
coefficients are getting multiplied by the sine or cosine of the concerned angle. This 
implies that in every succeeding steps, these coefficients are decreasing in their magni- 
tude. To justify the above statement, let us now consider all the coefficients obtained 
in the second iteration. 

At the end of stage 1 of the second iteration, the rotation coefficients become 


2f 21 = 2f 13 cosa 2 = ^gusinPjSinyiCOSOtj 
2g 21 = -2f 13 sina 2 = 2g u sinp 1 siny 1 sina 2 

bi 3 - 3 j 3 

2h 21 = 0 where cot2a 2 h 13 

At the end of the second stage of the second iteration, the rotation coefficients become 

2 {22 = 0 where cot2p 2 = — — f 21 

2 g 22 = 2g u sinp 1 sirry 1 sina 2 cosp 2 

Zh^ = -2g n sinp 1 sinyisina 2 sinp 2 

Similarly at the end of the final stage of the second iteration, the rotation coefficients 
reduce to 


2f 23 = -2g[ j sinP^in^sina^inPjsinyz 
b] 3 — 3 j3 

2 g 2 3 = 0 where cot2a 2 = h 13 

2h 23 = -2gj 1 sinp 1 siny! sina 2 sinp 2 cosy 2 

a 2 , p 2) and y 2 are the respective rotation angles along the z, x, and y axes in the 
second iteration. Hence it is observed with each iteration that the rotation coefficients 
get smaller and smaller in magnitude and eventually drop out. 

We are now in a position to formulate a rotation matrix whose elements 
correspond to the directional cosines of the x, y, and z axes of the rotated object. 


Rotation Matrix = Rp R a 


( 19 ) 


where 


R« = 


cosa sina 0 
-sina cosa 0 
0 0 1 


1 0 0 


Rp = 


0 cosP sinP 
0 -sinP cosP 


cosy 0 siny 


R r = 


0 1 0 
-siny 0 cosy 


Subsequently, 


where 




cosacosy + sinasinPsiny 
-cosPsina 

-siny:osa + cosysinasinp 


cosysina - sinysinpcosa sinycosp 
cospcosa cosP 

-sinasiny-cosysinPcosa cosPcosy 


( 20 ) 


a* = to,, p* = Xp b and y* = . 

i=l i=l i=l 

n corresponds to the iteration where all the rotation terms go to zero. 

Once the rotation terms, i.e., xy, yz, and xz are eliminated, the 3-D surface has 
the representation of 


F(x,y,z) = Ax 2 + By 2 + Cz 2 + 2Px + 2Qy + 2Rz + D = 0 (21) 

where A, B, C, P, Q, and R are the coefficients evolved after the elimination of the 
rotation terms. 

A natural question to ask is the following : " Can the terms of the first degree be 


proved numerically and experimentally. 


eliminated by means of a translation ?"The answer is, "sometimes they can and some- 
times they cannot". For the case, where the term can be eliminated, is supported by the 
following theorem [11]. 

3.2. Translation of the rotated object 

Theorem: The terms of the first degree of an equation of a quadric surface can be 
eliminated by means of a translation if and only if the surface has a center, in which 
case the first degree terms are eliminated if and only if the new origin is a centerfll]. 

Since all the 3-D object we are considering do have centers, we do not have to 
worry about the second case. The method of completing squares is the easiest to 
determine the coordinates of the new origin. 

Consider equation (21), grouping the like terms 


Ax 2 +2Px +By 2 + 2Qy + Cz 2 + 2Rz + D = 0 => 


F % 


f % 


r > 

X 2 + 2P-7- 
A 

^ J 

+ B 

y2 + 2Q b" 

+ c 

>. j 


Upon completing squares, we get 


P 

2 

r y 

Q 

2 

r •y 

R 

2 

P 2 Q 2 R 2 

X+ A 

b J 

+ B 

y + — 
J B 

b j 

+ c 

z+ c 

L J 

+ D - 

A + B + C 

< J 


where -P/A, -Q/B, and -R/C are the coordinates of the new origin. 

4. CONCLUSIONS 


( 22 ) 


In this research we have proposed a method to determine the pose and orientation 
of a natural quadric surface from its range image and subsequently eliminate all of 
these rotation parameters. 


Once the rotation and the translation procedures accomplish their objective of 
placing a quadric surface in a desired coordinate system, we will be in a situation to 
implement the necessary recognition process. We wish to extend upon our previous 
work [7] to the various quadric surfaces as mentioned before. Experiments for the 
alignment algorithm were performed on range data of circular cylinder rotated in 
space. The results obtained were very promising. Range data of quadric cones, para- 
bolic cylinders, and hyperboloids are presently being investigated. 
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